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Abstract 

We examine the adiabatic approximation in the study of a relativistic two-body problem with the gravitational 
radiation reaction. We recently pointed out that the usual metric perturbation scheme using a perturbation of 
the stress-energy tensor may not be appropriate for study of the dissipative dynamics of the bodies due to the 
radiation reaction. Over a time scale during which the usual perturbation scheme is valid, the orbits may not deviate 
substantially relative to the orbits of the background orbits. As a result, one can eliminate the orbital deviation 

*vO , through a gauge transformation. This is called the gauge problem of the gravitational radiation reaction exerted on 
the bodies, and it has been reported that a careful gauge fixing may be necessary to produce a physically reasonable 
prediction for the evolution of the system. 

We recently proposed a possible approach to solve this problem with a linear black hole perturbation. This paper 
proposes a non-linear generalization of that method for a general application of this problem. We show that, under a 

, & • specific gauge condition, the method actually allows us to avoid the gauge problem. 
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I. INTRODUCTION 



> 

^\ ' Many-body problems are fundamental problems in general relativity and have a long history of theoretical investi- 
t-H , gation. Starting from the famous paper of Einstein, Infeld and Hoffman, [1] the equation of motion of point particles 
■ was studied by various authors. The pioneering work of Einstein, Infeld and Hofman [1] used Dirac's delta functions 
as particles (the point particle approximation), and assumed a coordinate system in which the spacetime is weakly 
curved and the motion of the particles is sufficiently smaller than the light speed (the post-Newtonian approximation). 
To leading order in this approximation, we obtain a Newtonian equation of motion. Aside from the point particle 
"»«. . approximation, one can perturbatively consider higher-order terms of the post-Newtonian approximation, and Ein- 
stein, Infeld and Hofman derived the equation of motion with the leading relativistic effect under the post-Newtonian 
I , approximation. 

Much effort has been dedicated to deriving an equation of motion for a two-body system in a higher order of this 
approximation. [2] Such attempts were motivated by recent ground-based gravitational-wave interferometric detection 
projects, such as LIGO, VIRGO, TAMA and GEO. [5] These detectors are expected to facilitate observation of 
k> | gravitational waves from a close binary system consisting of stellar mass compact objects and, in order to extract 
astrophysical information concerning the gravitational-wave source from the observed waves, it is necessary to possess 
a theoretical gravitational waveform from a binary system. Such a binary is called an equal mass binary system and 
is believed to be accurately approximated by a relativistic two-body system, because the dynamical effect caused by 
the environment may not be important. Using the gravitational- wave generation formalism, [2] one can calculate the 
gravitational waveform from a two-body system, whose dynamics are governed by the relativistic equation of motion. 
We expect that detected gravitational waves might reveal a fundamental property of the gravitational law in the 
strongly gravitating limit. [5] 

The study of the gravitational waveform for a binary has also been carried out using black hole perturbations. [6] In 
this case, we consider the case that the mass ratio of the binary is extreme, such as in the case of a binary consisting 
of a supermassive black hole and an inspiralling stellar mass compact object. It is expected that gravitational waves 
from such a binary system will be detected by a space-based gravitational-wave interferometric detector, such as 
LISA. [5] Because of the extreme mass ratio, the gravitational potential of the system is due almost entirely to the 
heavier body, and we may suppose that the heavier body creates a background geometry, which is assumed to be a 
black hole geometry. The gravitational effect caused by the lighter body is treated as a metric perturbation of the 
background black hole geometry. Given the motion of the lighter body, the black hole perturbation formalism yields 
the gravitational waveform from the binary system, [6] with which we may study the gravitational theory. 
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During the observation time in both cases, the gravitational radiation reaction becomes important for the orbital 
motion of the binary. In this paper, we discuss the problem of the mathematical framework for the study of this 
effect. Until this time, in most cases, the gravitational radiation reaction effect has been dealt under the so-called 
adiabatic approximation. There are two key assumptions for the adiabatic approximation. The first approximation 
is that the instantaneous dynamics of the orbit are conservative, because the radiation reaction effect is weak. We 
specifically consider the case in which an orbit is characterized by several constants of motion, denoted by £ a , and 
we consider the evolution of such "constants" of motion instead of the orbital coordinates of the bodies. With this 
approximation, we calculate the radiation reaction exerted on the "constants" of motion by the conservative orbit of 
given constants of motion. Thus, we have an evolution equation of the form 

(1.1) 

The second assumption is that the time averaged change of the "constants" of motion < d£ a /dt > describes the 
orbital evolution, because the radiation reaction affects the orbital dynamics in a secular manner. The "constants" 
of motion evolve in the so-called radiation reaction time (T rac j), and it is believed that the adiabatic approximation 
holds when the radiation reaction time is much longer than the dynamical time scale of the orbit (Tjy n ). 

In the case of an equal mass binary, we consider a nonlinear metric perturbation with a flat metric as the background. 
[2] We separate the metric perturbation into the symmetric and anti-symmetric parts with respect to the background 
coordinate time. The conservative motion of the bodies is derived using the symmetric part of the metric, and it is 
derived simultaneously and consistently with the symmetric part of the metric under the post-Newtonian expansion. 
[2] When the spins of the bodies are negligible, the conservative motion of the binary is characterized by two constants 
of motion, the orbital energy and the angular momentum. Using the conservative motion of the bodies as a source, one 
can derive the anti-symmetric part of the metric with the wave generation formalism, employing the post-Minkowskian 
expansion. [2] This part is responsible for the gravitational radiation reaction and, by evaluating the gauge-invariant 
flux of the gravitational radiation momentum escaping to infinity, [3] we obtain the time averaged losses of the orbital 
energy and the angular momentum as functions of the constants of motion. 

In the case of the extreme mass ratio binary, we take the evolution of a geodesic of the background black hole 
geometry as a conservative motion of the orbit. When the background is a Kerr black hole, the geodesic is characterized 
by three "constants" of motion, the orbital energy E, the z-component of the angular momentum L z , and the Carter 
constant C. The metric perturbation induced by the geodesic is responsible for the gravitational radiation reaction. 
The radiation reaction exerted on these "constants" is derived by evaluating the radiative part of the self-force, * and 
we derive a formula for the time averaged losses of these three "constants" as functions of the constants of motion. It 
is also shown that these time averaged losses are gauge invariant in Ref. [4] . 

Thus, we have successfully formulated the radiation reaction problem in both cases under the adiabatic approxima- 
tion, and the results in both cases are shown to be consistent in the post-Newtonian limit and in the extreme mass 
ratio limit. [2,6] 

The adiabatic approximation is merely an approximation for a practical calculation of the gravitational radiation 
reaction effect, and it would eventually be invalid for an accurate prediction of gravitational waveforms. In this 
paper, we examine the validity of the adiabatic approximation by estimating the effect on the orbital evolution which 
appears when the adiabatic approximation becomes invalid. Regarding the second assumption, an explicit calculation 
of the part of the self-force other than the time averaged losses of the "constants" was done in the equal mass binary 
case. [9] In the extreme mass ratio binary case, we showed that the time averaged losses of the "constants" are the 
only gauge invariant components of the self-force, and that one may take the gauge transformation to eliminate the 
gauge dependent components on the radiation reaction time scale. In this paper, we study the first assumption of the 
adiabatic approximation. To this point, we have discussed the adiabatic approximation of the orbit. In Ref. [7], we 
pointed out that for a consistent study of the validity of the adiabatic approximation, it is necessary to consider the 
adiabatic expansion of the metric together with the adiabatic evolution of the orbit. We carried out a preliminary 
investigation of the adiabatic expansion of the metric in Ref. [7] and here we extend the idea proposed there. 

The organization of this paper is as follows. In §11, we briefly summarize our present result concerning the validity 
of the adiabatic approximation. We have previously carried out several investigations of the validity of the second 



^It is believed that the self-force has a conservative component. However, in Ref. [7], we find that one can eliminate this part 
by taking the radiation reaction gauge condition over the radiation reaction time. This shows that the conservative component 
of the self-force is not observable unless we can observe gravitational waveforms with a precision of the mass ratio in each 
period. 
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assumption. However, the validity of the first assumption has not been carefully studied until now. In §111, we argue 
that the question concerning the first assumption is related to the basic mathematical scheme used to treat a binary 
system coupled to a gravitational field. In §IV and V, we propose a new metric perturbation scheme and its practical 
calculational procedure appropriate for studying the evolution of the metric and the orbit. In §VI, we show that the 
resulting perturbation expansion of the metric behaves properly under a specific gauge condition on the second-order 
metric perturbation. We summarize results in §VII. 

II. ADIABATIC APPROXIMATION 

As pointed out in the previous section, the adiabatic approximation uses two separate assumptions. Here we 
summarize our present result on the validity of these two assumptions separately. 

The first assumption is to use the conservative orbit as a source to calculate the radiative component of the 
gravitational field. The conservative orbit is believed to be a good approximation of the orbit instantaneously, but, 
the orbit will deviate from a conservative orbit due to the radiation reaction. In order to examine its validity, it 
is necessary to consider the component of the gravitational field induced by the orbit which evolves due to the 
radiation reaction. At this point, it is important to note that the orbital coordinates without the radiation reaction 
are characterized by not only the constants of motion, £, but also the so-called phase constants, C (see §111). While 
the effect of the evolving "constants" of motion on the metric perturbation has been studied, the effect of the evolving 
phase "constants" had not been carefully considered until we pointed it out in Ref. [7] . This is because one can set the 
phase constants to zero for a conserved orbit through an appropriate isometry transformation, but once we consider 
the orbital evolution due to the radiation reaction, one cannot set the phase "constants" to zero during the evolution. 
Because the evolution time scale for the phase "constants" is shorter than that of the "constants" of motion, it is 
more important to consider the evolution of the phase "constants" . [7] 

A consistent calculation of such a metric can be done by calculating the evolution of the metric and the orbit 
together. In Ref. [7], we propose an adiabatic approximation for both the metric and the orbit in the case of an 
extreme mass ratio binary system, and we find that one can predict the orbit under this first assumption over the 
radiation reaction time scale. 

Because the radiation reaction time scale is approximately equal to or less than the observation time scale in the case 
of an equal mass binary for an observation by ground-based detectors and in the case of an extreme mass ratio binary 
for an observation by a space-based detector, the result given in Ref. [7] may pose a serious question concerning the 
reliability of theoretical waveforms derived using the adiabatic approximation. In §111, we discuss the consequences 
of this problem in greater detail, and here we merely comment that, while we have been considering a method to 
construct a metric accurately at an instance during the evolution due to the radiation reaction under the adiabatic 
approximation of the orbit, we focus on how to construct an accurate evolving metric together with the orbit. An 
accurate derivation of the metric and orbit at an instant is a necessary step, but this does not necessarily give us the 
accurate evolving metric and orbit. We refer to this as the "Burke problem". 1 '' 

The validity of the second assumption was recently investigated in Ref. [9] for the case of an equal mass binary, 
and in Ref. [7] for the case of an extreme mass ratio binary. The correction to the second assumption is called the 
"phasing effect" . The evolution equations of the "constants" of motion, (1.1), become generally time varying functions 
for fixed £ a . With the second assumption of the adiabatic approximation, we usually take only the part of (1.1) that 
does not depend explicitly on time. Thus, a correction to the second assumption can be obtained by calculating the 
full evolution equation (1.1), which depends explicitly on time. 

This problem was not considered until recently in the case of an equal mass binary system, because this effect does 
not appear when the relative orbit of the two body system is instantaneously circular and because the orbit of an 
equal mass binary system is circular to a very high precision in the observationally most interesting case. [10] When 
the orbit is circular, the metric and the matter source allow the helical Killing vector, £ a = (d/dt) a + Q(d/d(j)) a 7 
where f2 is the angular velocity of the orbit for a broad range of choices of gauge conditions. The 4- velocity of the 



^Professor William Burke was a pioneer of the gravitational radiation reaction problem. When he was a student at CalTech, 
he challenged his supervisor, Professor Kip Thorne, with this question. One can see a record of Burke's challenge in front of 
Thome's office at CalTech. Although he conceded by his famous discovery of the so-called Burke- Thorne radiation reaction 
potential [8] under the post-Newtonian approximation, the question of mathematical validity involved in the treatment of the 
radiation reaction problem remains a fundamental question in the study of the gravitational radiation reaction. 
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orbit is proportional to the helical Killing vector, and the radiation reaction acting on the orbit is constant along this 
Killing vector. Thus, (1.1) does not depend explicitly on time for a circular orbit. 

This problem has an aspect to which a special attention must be paid. Although the time averaged components 
of (1.1) are gauge invariant, the rest of (1.1) is gauge dependent in general. § Thus, we need the additional argument 
that those gauge dependent components are actually what we can observe by calculating gravitational waveforms at 
the spatial infinity. 

To conclude this section, we summarize the qualitative behavior of the orbital evolution due to the radiation 
reaction without making the second assumption of the adiabatic approximation. [7] We suppose that the adiabatic 
metric perturbation in Ref. [7] can be defined in the equal mass binary case and that the result of the radiation 
reaction induced by a "conservative orbit" can be used to infer the evolution of the "constants" of motion over the 
radiation reaction time scale. The "constants" of motion evolve qualitatively as 

£ a = £ a Jl + -^ + 5] , (2.1) 



-^rad 

where £§ are the initial values, and T ra j is the radiation reaction time defined by T rac | = £/ < d£/dt >.** The 
second term on the RHS of (2.1) comes from the time averaged components of (1.1), and the third term comes from 
the rest of (1.1). If we assume that the amplitude of (1.1) is 1/T rac j and that (1.1) varies on the dynamical time scale 
of the orbit, 7jy n , we obtain the relation 5 ~ ^dyn/^rad" 

This is not the end of the story, however, because what we observe is the phase of the gravitational waves, which is 
twice the angular position of the binary. Noting that the instantaneous frequency, fi, is a function of the "constants" 
of motion, the phase evolves as 

$(t) = / dm(£ a ) = n(£«) ( t + -J— + s p ) , (2.2) 

J V J p.rad / 

where rac j = (2fl/ d a fl£Q )7 1 ra( j ~ T ra( ^. Because of the time integration, we have S p ~ ^dyn^/^rad" 

Equation (2.2) suggests that S p becomes important at the radiation reaction time scale, t <~ T p ra j for an accurate 
prediction of the phase, and we need to consider the phasing effect, which is ignored when the second assumption of 
the adiabatic approximation is employed. In the case of an equal mass binary, it is necessary to calculate the phasing 
effect with an appropriate gauge condition to obtain an accurate prediction of the gravitational wave phase over the 
radiation reaction time scale. However, as we discuss in Ref. [7], in the case of an extreme mass ratio binary, we may 
choose a gauge condition under which we have 5 P ~ /T ig £, thus, the phasing effect is not important for a reliable 

prediction of gravitational waveforms. 



III. GAUGE PROBLEM AND FIRST PRINCIPLES CALCULATION 



With the first assumption of the adiabatic approximation, we consider first calculating the instantaneous metric 
induced by a conserved orbit as accurately as possible, and then deriving the evolution equation for the "constants" 
of the form (1.1) from the instantaneous metric. In order to ascertain the validity of the first assumption, one has to 
estimate the extra terms coming from the evolving metric induced by the orbit deviating from a conservative orbit 
due to the radiation reaction. We refer to the derivation of (1.1) together with a consistent derivation of the metric 
as the "first principles calculation" . 

The most important problem here is properly choosing the gauge condition. The gauge condition entirely determines 
the orbital evolution, and it has been shown that one may even choose a gauge condition such that the RHS of (1.1) 
vanishes consistently with an usual metric perturbation scheme. [7] However, such a choice of the gauge condition 
is not appropriate, because the perturbative calculation of the metric becomes invalid on a short time scale, and 



^Strictly speaking, even the time averaged components of (1.1) are gauge dependent. However, in a sufficiently large and 
physically reasonable class of gauge conditions, we showed in Ref. [7] that the time averaged components of (1.1) become 
invariant. 

"Precisely speaking, the expressions given in (2.1) and (2.2) are valid only when t <ti T mc j. Terms which behave as ~ t 2 /T^ ^ 
appear in (2.1) when t ~ T rac j. However, it has been argued that as long as we are interested in the validity of the adiabatic 
approximation, its qualitative features can be understood sufficiently well from (2.1) and (2.2). [7] 
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consequently, one cannot make any meaningful prediction of gravitational waveforms modulated by the radiation 
reaction. [7] Unfortunately, we have not yet found a concise criterion for this gauge condition. In order to determine 
whether or not a gauge condition is appropriate to study the radiation reaction problem, one must calculate the 
evolution of the metric together with the orbital evolution due to the radiation reaction. If the resulting evolution 
of the metric is an accurate approximation for a long time (longer than the radiation reaction time scale), we can 
conclude that the gauge condition is appropriate to study the radiation reaction problem on this time scale. 

Well known methods for solving the Einstein equation, such as the post-Newtonian approximation for an equal 
mass binary and the black hole perturbation for an extreme mass ratio binary, are applicable when the binary motion 
is conservative. Therefore, one can derive the instantaneous metric with these methods. The evolving metric may 
be derived by applying the adiabatic expansion to this instantaneous metric, as we proposed in Ref. [7]. Although 
the original idea of the adiabatic metric perturbation given in Ref. [7] seems successful for the study of the radiation 
reaction of an extreme mass ratio binary, it was defined only for a linear perturbation, and thus, it is not clear whether 
the idea can be extended to the case of an equal mass binary for which the non-linear calculation of the metric with a 
flat background is done under the post-Newtonian approximation. In addition, the equation of motion resulting from 
the adiabatic metric perturbation poses a serious question with regard to its predictability. 

Using the adiabatic approximation, we can calculate the instantaneous metric induced by a conserved orbit with a 
certain approximation, such as the 2nd post-Newtonian approximation in the case of an equal mass binary and the 
black hole linear perturbation in the case of an extreme mass ratio binary. The orbital equation for the radiation 
reaction (1.1) is derived from the instantaneous metric. Recall that we restrict our consideration to the case that the 
4-velocity is determined by the constants of motion, £. The orbital coordinates as functions of the orbital parameter 
can be obtained by further integrating the 4-velocity, and the integral constants for this integration are called the 
"phase constants", C. Thus, the conserved orbit is determined by £ and C, and the instantaneous metric is derived as 
a function of £ and C. As a result, the RHS of (1.1) becomes a function of £ and C under the adiabatic approximation 
in general. Applying an appropriate isometry transformation, one can make the phase constants vanish in the case of 
an equal mass binary. We showed that this is true for the extreme mass ratio binary case in Ref. [7]. Therefore, under 
the adiabatic approximation of the orbit, the RHS of (1.1) is a function of only £, and one can derive the evolution 
of £ as a function of the orbital parameter in a closed form using (1.1). 

Now, let us consider a correction to the adiabatic approximation. We consider the orbital evolution due to the 
radiation reaction; that is, £ and C are no longer constant. Then, it becomes necessary to consider the explicit 
dependence on C of the metric, because the isometry transformation cannot eliminate C during the evolution. Fur- 
thermore, the dependence on C is more important than that on £, because C evolves more rapidly than £ . [7] In Ref. 
[7] , preliminary discussion on this effect is given for the case of an extreme mass ratio binary. We define the adiabatic 
linear metric perturbation as that consisting of the replacement of the constants £ and C by the evolving £ and C 
of the instantaneous metric. Through this operation, the evolution equation (1.1) becomes that which we use under 
the adiabatic approximation of the orbit, to leading order in the perturbation expansion, and it does not depend on 
C. Next we consider the equation for the correction to the adiabatic linear metric perturbation. Although we could 
not obtain a consistent expansion scheme in Ref. [7] , the source term for the correction to the adiabatic linear metric 
perturbation includes terms proportional to the first and second derivatives of £ and C. If we consider the linear 
correction of the adiabatic linear metric perturbation, it may be reasonable that the correction metric includes terms 
proportional to the first and the second derivatives of £ and C. In this case, instead of (1.1), the orbital equation 
becomes 

rlfa _ 

— =F a (£,£,£,C,C), (3.1) 

and one cannot solve the orbital evolution for £ in a closed form. 

One can eliminate the dependences on £, £, C and C using the orbital equation derived with the adiabatic linear 
metric perturbation. In Ref. [7], we derived the qualitative behavior of the derivatives of £ and C employing the 
adiabatic linear metric perturbation as 

£ ~ £^J—, £ ' 1 



T rad T rad 

C ~ — (E — En) , C r * J — E , 

d£ { 0j ' d£ T rad ' 

where is the frequency associated with the phase constants, which is a function of £, and £ is the initial value of 
£ when we switch on the radiation reaction. Because the behavior of C involves £ , we have two serious problems, 
described below. 
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1) Adiabaticity: The orbital prediction obtained this equation depends on the choice of the initial value, £q. 

2) The Burke problem: Because C grows in time, the correction to the adiabatic linear metric perturbation may also 
grow in time. As a result, the correction would eventually become dominant over the adiabatic metric. 

The hrst problem suggests that we might not be able to predict the gravitational waveform, because the result depends 
on the choice of the initial time. This would imply that our basic hypothesis of the adiabatic evolution of the orbit 
might be wrong. The second problem suggests that the adiabatic linear metric perturbation may not be a good 
approximation of the metric at large times. 

These problems are unexpected with our current understanding of the gravitational radiation reaction problem. 
Actually, it is widely accepted by the physics community that the orbital evolution calculated under the adiabatic 
approximation was confirmed by a Nobel Prize winning observation of the Hulse- Taylor binary pulsar. [11] 

We believe that the cause of these problems is an inappropriate choice of the gauge condition. As we stressed in Rcf. 
[7], the calculation could become highly nonperturbative with an inappropriate choice of the gauge condition when 
we consider the evolution of the metric and the orbit due to the radiation reaction. For this reason, we cannot predict 
the gravitational waveform for a sufficiently long time. This is the Burke problem. If the Burke problem arises, the 
prediction is valid for only a short time, and the initial time can be chosen only within this short time span. This is 
the adiabaticity problem. We note that the source term for the correction to the adiabatic linear metric perturbation 
depends on the gauge condition. We hypothesize that an appropriate gauge condition is the condition that we can 
eliminate the C-dependence of the source term and we may be able to avoid these problems. 

In the following two sections, we propose a consistent perturbation scheme to construct the adiabatic metric in 
a nonlinear manner by extending the idea proposed in Ref. [7] This scheme consists of two different perturbation 
methods with the same small expansion parameter A. In §VI, we show that, under a specific gauge condition, one can 
actually avoid the Burke problem for the leading-order correction to the adiabatic linear metric perturbation. 

IV. RENORMALIZED METRIC PERTURBATION 

We consider the Einstein equation, employing the geometrical units, in whichG = c = 1. 

G a ^ = 8irT a P . (4.1) 

In this section, we discuss the formal application of the renormalized perturbation method to a metric perturbation, 
and we do not specify the form of the stress-energy tensor. We consider a general vacuum metric g( bfc ) Q ' 3 as the 
background for the perturbative expansion, so that one can use this scheme for both the equal mass binary case and 
the extreme mass ratio binary case. Below, the covariant derivative is taken with respect to the background. 

The expansion parameter A is taken to be the amplitude of the stress-energy tensor, T"' 3 ~ 0(A). Because the 
metric perturbation induced by this stress-energy tensor creates a radiation reaction that influences the evolution of 
the matter, we can assume T ali ' .p ~ 0(A 2 ). 

We expand the metric in A as 

The Einstein tensor can be schematically expanded as 
= A [G a ?] (1) (/i«) 

+A 2 | rga/J] (1) (/l (2) } + (2) (fcW 

+A 3 { [G Q/3 ] (1) (h^) + 2 [G a/3 ] (2) (h™, ftW) + [G Q/3 ] (3) (hW,hW,hW)} 

+ •••, (4.3) 

where [G af3 ] {1) (h) is linear in h at \ [G af3 ] (2) (h^\h^) is bi-lincar in both h^"! 3 and h^"! 3 \ and so on. We recall 

that there is an algebraic relation [G Q/3 ] ^ (h) = for any h a ^ . 

To obtain a consistent perturbation expansion, we could consider the expansion of the stress-energy tensor, T Q/3 , 
in A. If we do this, the stress-energy tensor is conserved in the background to leading order due to this consistency, 
and the dissipative dynamics of the matter due to the radiation reaction can only be taken into account at higher 
order. Because the dissipative effect accumulates in time, the higher-order terms eventually become dominant over 
the leading term, and the expansion of the stress-energy tensor soon becomes invalid. This implies that the Burke 
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problem is inevitable, and we cannot make a physically meaningful prediction of the system's evolution due to the 
radiation reaction. [7] For this reason, we do not expand the stress-energy tensor and propose a new perturbation 
scheme with which the dissipative dynamics of the stress-energy tensor can be renormalized in the leading order of 
the metric perturbation in a consistent manner. For this purpose, we use a gauge fixing term as a counterterm for 
the renormalized perturbation. This idea is similar to that of the so-called fast motion approximation. [12] However, 
the key difference is that in the present case, we have order-by-order equations for the metric perturbations at each 
order in the adiabatic expansion we propose in the next section. 

It is important to note that the stress-energy tensor does not satisfy the conservation equation in the background, 
T a/3 -/3 ^ 0, due to the radiation reaction caused by the metric perturbation. Therefore, it is inconsistent to simply 

equate A [G af3 ] (1) (ftW) with T a(i in order to define h^ al3 . We also note that, for a given divergence-free tensor 

ya/3 (f a P.p — o) ; the linear Einstein equation [G"* 3 ] ^' (h) = f a @ with an appropriate boundary condition does not 
have a unique solution, because of the gauge freedom. This is equivalent to the fact that the linearized Einstein 
operator [G"' 3 ]^ docs not have an inverse (the Green function) until we add a gauge fixing operator. We formally 

define the gauge fixing operator A a/3 7 s, which consists of differential operators that are at most second order. We 
have a tensor Green function g a/3 -y'S' {x, x') that satisfies 

{[G^ + [A"%}. 9 ^,(x,x') 

= V^^^) 6 ^^^ (4-4) 

w here x — x 1 ) is the 4-dimensional Dirac delta function and </ bfe ) is the matrix determinant of g„p* ■ We choose 
the gauge fixing term such that the linearized Einstein equation is no longer divergence free in an algebraic manner 
and, as a result, it allows a general symmetric tensor as a source. We remark that the Green function defined in this 
manner includes propagation of not only the gravitational mode but also the vector and the scalar modes. 
We define the linear metric perturbation Xh^ al3 of (4.2) as 

Xh (1 ^ = J d 4 x'^-g( bk )g a P Y5 > {x, x')SnT c ' d ' (x') . (4.5) 
By definition, the linear metric perturbation satisfies 

[G aP ] (1) (Xh {1) ) = 8ttT q/3 - A Q/3 (Ah (1) ) . (4.6) 

We recall that T af3 '.p ~ 0(A 2 ) and we require that Xh^ af} be a solution of the linear Einstein equation to 0(A). 
Then, we have the condition 

A a/3 (A/i (1) )~0(A 2 ). (4.7) 

Taking the divergence of (4.6), we have 

T"% = J-A^A/iW) ~ 0(A 2 ) , (4.8) 

which gives the equation of motion to an accuracy of O(A) as a consistency condition. 

The gauge fixing term of (4.6) must be taken into account as a source for the second-order metric perturbation, 
together with the quadratic contribution of the linear metric perturbation. Using the tensor Green function, we define 
the second-order metric perturbation as 

X 2 h^ af3 = J d 4 x'g af3 7l5 >(x,x') 

x (^' 5 ' (Xh^){x') - [G 7 ' y ] (2) (Xh^,Xh^)(x'^j . (4.9) 

By definition, the second-order metric perturbation satisfies 

[G^] (1) (A 2 ^ 2 >) + [G^] (2) (XhW,XhV) 

= K af) (Xh { ^)~K af \X 2 h^). (4.10) 
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As (4.7), in order to make X 2 h^ al3 an approximate second-order metric perturbation to 0(X 2 ), we impose the 
condition 

A«/3( A V 2 >) ~0(A 3 ). (4.11) 

Taking the divergence of (4.10), we have 

T«% = ^A«%(A/,«) 

= _L [G^l (2) (Xh^,Xh^). fj + 0(X 3 ) . (4.12) 

This shows that the equation of motion to 0(A 2 ) can be derived with the linear metric perturbation \h^ al3 . 
In a similar manner, we define the n-th order metric perturbations as 



x(A 7 ' 6 '(A/i ( "- 1) )(a;')-2 [G 7 ' 5 'l (2) (A"" 1 /^"" 1 ), A/i (1) )(a/) 



(n) 



(A/iW- ••)(*')) • (4-13) 



This satisfies the equation 



[G a/3 ] (1) (A n ftW) + • • • + [G a/3 ] (n) (A/i^ • • •) 

= A afl {X n - 1 h {n - 1) )- K af3 {X n h {n) ). (4.14) 



By imposing the condition 



A«/3(A"/t(")) ~0(A" +1 ), (4.15) 
we obtain the equation of motion to 0(A") as 

+2 [G^} (2) (X 2 h^,Xh^), fj + [G ab ] i3) (XhW,XhV,XhV). b 
+ ■■■ 

+2 [G ab ] (2) (A"" V"" 1 ), A/i«) ;6 + • • • + [G ah ] (n) (A/*« , A/i« , • • -) ; & 
+0(X n+1 )] . (4.16) 



The equation of motion to 0(A") can be derived by calculating the metric perturbations, Xh^ al3 (k = 1, • • • , n — 1). 

We note that the gauge condition is obscure in this perturbation scheme. This is because we do not expand 
the stress-energy tensor in A, and as a result, the stress-energy tensor does not satisfy the conservation law in the 

background metric. We do not equate A [G"' 3 ]^ (h^ ) and T Q/3 , as in the usual metric perturbation scheme, and, as 
we see from (4.6), we incorporate the difference of these terms in the gauge fixing term. Therefore, the gauge fixing 
term does not vanish in this perturbation scheme. 

We also note that the validity of this perturbation scheme depends on the gauge fixing term. In order to see whether 
or not we can avoid the Burke problem considered in §111, we must explicitly check the behavior of the gauge fixing 
term. In §VI, we consider this problem with a specific gauge fixing term. 

V. ADIABATIC EXPANSION 

Following the formal discussion of the metric perturbation scheme given in the previous section, we define the metric 
perturbation which includes the radiation reaction effect to leading order in a consistent manner as (4.5), (4.9) and 
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(4.13). We may still have a technical difficulty arising in an explicit calculation of these metric perturbations. Because 
the source term for h^ a ^ cannot be derived until we derive the equation of motion to a necessary order, the metric 
perturbation scheme presented in the previous section is not systematic in the sense that one cannot derive a metric 
perturbation of a given order from lower order metric perturbations. A similar difficulty appears in the so-called 
fast motion approximation. [12] In order to avoid this technical problem, we further use another perturbation scheme 
proposed in Ref. [7] which uses the same small parameter, A. For simplicity, we only discuss the derivation of h^ al3 
by integrating (4.6). We consider the application of the adiabatic expansion proposed in Ref. [7] with h^ af3 of the 
form 

Xh^ = Xh^ afi + AV 1 ' 1 ^ + • • • . (5.1) 

We suppose that we are able to calculate the higher-order metric perturbations (4.9) and (4.13) in the similar manner 
through the expansion 

X n h (n) a f3 = x n h (nM)afJ + A n+l ft (n,l)a/3 + . . . _ (5 2) 

We consider the stress-energy tensor of point particles and restrict ourselves to the case in which the equation of 
motion (4.16) is integrable, with constants of motion £ a in the absence of a gravitational radiation reaction, as in the 
case of a binary system of non-spinning particles. As discussed in §111, we also need the phase constants C b in order 
to specify the orbital coordinates, and the conserved stress-energy tensor is described as a function of £ a and C b . We 
suppose that there is a method to calculate (4.5) for a conserved stress-energy tensor. Using this method, we have the 
"ordinary" linear metric perturbation h^ a ^(x;£ a ,C b ), which is the solution of the usual linearized Einstein equation 

A [G ati ] {1) (U iy ) = T a P(x;£ a ,C b ). (Because we do not consider the radiation reaction at this stage, the stress-energy 
tensor on the RHS is conserved in the background.) 

Now we consider the effect of the radiation reaction. Assuming the adiabatic evolution of the orbits, we consider the 
"constants" as functions of the local time t which evolve due to the radiation reaction. Assuming the small parameter 
A used in §IV, we have ^£ a ~ ^C b ~ 0(A). In Ref. [7], we extend this picture of the adiabatic evolution of the orbit 
to the evolution of the linear metric perturbation. We foliate the spacetime into smooth spacelike hypersurfaces, and 
the foliation function, f(x), is taken such that it coincides the local time t around the particles. The adiabatic linear 
metric perturbation is defined by using the "ordinary" linear metric perturbation of the constants, £ a (f) and C b (f), 
on the foliation surface, f(x) = f. Because liniA^o (T^f / A) — > 0, the adiabatic metric perturbation gives (4.5) to 
0(A). Therefore, we take the adiabatic linear metric perturbation as the leading term in the adiabatic expansion of 
(5.1): 

h (i,o )a p = ~ h W"0( X ; £ a (f(x)), C b (f(x))) . (5.3) 
This yields the orbital equation (4.12) to 0(A) as 



dt 



—£ a 
dt 



d 



— C b 

dt 



(1) +0(A 2 ). (5.4) 



The equation for the correction to (5.3) can be obtained by inserting (5.1) into (4.6). Because (4.6) is linear in 
/j(i)"/3 we have 



{ [G aP ] (1) (A n /i< 1,n >) + A Q/3 (A n /i (1,n) )} 



71=1,2,- 

= - [G aP ] {1) (Xh (h °'>)-A aP (Xh^^)+T a ^(x;£ a (t),C b (t)). (5.5) 

We recall that the linearized Einstein operator and the gauge fixing operator consist of differential operators that are 
at most second order. Therefore, the RHS of (5.5) becomes 

\ ( c upd£ a d£ a aj3 d£ a dC b af3 dC a dC b \ . . , . 

y a "'~dT~dT + s ^ b ^r^T + bb ~dT~dT ) v ^ v ^) 



( r l2ca rl 2 r b \ 
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where d n £ a /dt n and d n C b /dt n are evaluated at t = f(x). Using (5.4), we obtain the following equation for ft/ 1 ' 1 )"/ 3 
by taking the terms of 0(A 2 ) of (5.5): 



[G af3 ] {1) (aV^ha^aV 1 - 1 )) 



A \ Sa dt 



dt 



£" 



(i) 







(i) 










d c b 




dt 




dt 




(5.7) 



Thus, we have the equation for ft/ 1 ' 1 )"/ 3 , where source term is derived using the metric perturbation ft/ 1 ' )"/ 3 and the 
orbital equation of motion, [:§£ ]^ and [jr^ 6 ]^. However, (5.7) cannot be integrated, because ft/ 1 ' )"' 3 depends 
on £ a (t) and C b (t). We consider repeating the adiabatic calculation used for ft/ 1 ' )"* 3 in order to integrate (5.7). In 
other words, instead of solving (5.7) exactly, we integrate the equation, regarding £ a and C b on the RHS as constants. 
Then, we replace £ a and C b by £ a (f) and C b (f). By this method, we can solve (5.7) to an accuracy of 0(A 2 ), and we 
take into account the extra term of 0(A 3 ) coming from this adiabatic calculation of ft/ 1 ' 1 )"* 3 as an additional source 
term for ft/ 1 ' 2 )"/ 3 . 

The equation for ft/ 1 ' 2 )"/ 3 can be derived from the metric perturbations ft/ 1 - )"' 3 and ft/ 1 ' 1 )"/ 3 , and the orbital equa- 
tions [f t £ a ] (1 \ [£ t £ a ] {2 \ [ft Cb ] {1) and [li Cb ] {2) - In ordcr to derive [li £a ] {2) and [li Cb ] {2) ' ^ (1,1)Q/3 is insufficient, 
and it is necessary to calculate ft/ 2 ' ) Q/3 , as we see from (4.16) with n = 2. The equation for ft/ 2 - )"' 3 can be derived 
from ft/ 1 ' )"/ 3 . Thus, ft/ 1 ' 2 )"/ 3 can be derived from lower-order metric perturbations. 

Repeating this procedure, we find that the equation for ft/™'™)"/ 3 can be derived from the metric perturbations 
h (n',m')af3 <y + m ' < n + m ) an d the equations of motion [f t £ a ] ^ and [f f C b ] ( "° (n' < n + m). The equations of 

motion [jj£ a ]^™' and [^C h ]^™' > are derived from the metric perturbations ft(™/ m ')"/3 ( n ' _|_ m < < n y Thus, with this 
metric perturbation formalism, we can systematically calculate the metric perturbations order by order. 



VI. GAUGE CONDITION AND THE BURKE PROBLEM 



To this point, we have discussed the new metric perturbation scheme in §IV and §V without specifying the gauge 
condition, and there is a possibility that the Burke problem exists with an inappropriate gauge condition. In this 
case, although the expansion of the metric in the small parameter A is initially consistent, it would become invalid 
as the system evolves for a sufficiently long time after we switch on the radiation reaction. In this section, we choose 
a specific gauge fixing term in the new metric perturbation scheme and show that we can avoid the Burke problem 
with this gauge condition up to second order in the metric perturbation. 

We define the gauge fixing term as 

A af} 7 g = \ (Vv 53 f^ + V p V s g^ a - s (6i)a "V 7 vV) • (6.1) 
Then, the renormalizcd perturbation equations for ft( 1 ) Q/3 and ft( 2 ) Q/3 become 

^ ^(i)«/3;7. 7 + 2R a 7 P s h WlS ) = 8irT af) , (6.2) 

y (ft( 2 )«^ ;7 + 2i?VW 2 ^) 

= - [G«' 3 ] (2) (Aft( 1 ),Aft«) + A Q ' 3 (Aft( 1 )), (6.3) 

where the Riemann tensor, i?"-^, is defined for the background metric. Because (6.2) does not necessarily require the 
divergence-free condition for T a @ , the Green function of (6.2) includes the propagation of not only the gravitational 
mode but also the vector mode and the scalar mode as discussed in §IV. 

We first treat ft/ 1 ' 1 )"/ 3 by carrying out the adiabatic expansion of (6.2). Let us denote the Green function with 
an appropriate boundary condition by g a ^ 1 'S'(x,x'). We recall that the stress-energy tensor is a function of the 
"constants" evolving due to the radiation reaction, T Q/3 = T al3 (x; £(t),C(t)). In the adiabatic expansion presented in 
§V, ft( 1 '°) a ' 3 can be defined in terms of the Green function as 
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J dx'V^g^g^ys' (x, x')t^' s ' (x' ; £(f), c(/)) 



(6.4) 



/=/(*) 



Then, the equation for A 2 /i^ 1 ' 1 ) a/3 (x) becomes 
A 2 



^-(/ l ( 1 ' 1 )^ ;7 + 2i?V^ (1 ' 1)75 ) 

= 8 ttT^ - ^ (MW^ + 2i?V^ (W ) + 0(A 3 ) . 
Using (6.4), the RHS of (6.5) becomes 

J dx>^^g a P YS ,(x,x')±T'>' s '(x';£(f),C(f)) 



(6.5) 



-8tt 
-8tt 



/=/(*) 



| dx'y^Wg^s, (x, x')^' 5 ' (x';£(f), C 



(/)) 



/=/(*) 



(6.6) 



The first term has a part of 0(A 2 ) dependent on (d/df)C. However, by taking a smooth foliation, one can make 
f ;t - e (x) arbitrarily small, and we can eliminate this term. The second term also has a part dependent on {d/dj)C 1 
but it is of 0(A 3 ) and can be taken as the source term for /i( 1 ' 2 ) a/3 . Thus, the source term of (6.5) does not have an 
explicit dependence on the elapsed time of (d/df )C, and we avoid the Burke problem for /j( 1 ' 1 ) Q ' 3 . 

We next treat h^ 2,0 ^ ^ by performing the adiabatic expansion of (6.3). The equation for /j( 2 > ) Q ^ becomes 



= - [G Q/3 ] (2) (A^ 1 '°),A^ 1 '°)) + A Q ' 3 (A^ 1 ^) + 0(A 3 ). 



(6.7) 



The first term on the RHS of (6.7) has a term dependent on (d/df)C, but it is of 0(A 3 ). For the second term on the 
RHS, instead of using (6.4), we consider the equation for A a ^(\h^). By a simple manipulation, we obtain 



[A Q/3; ^ ;7 + 2R a ^ 5 ^ S ] (A/i (1) ) 



= - 8tt T a T.~^ + 



rpflj :a _ 1 (6fe)a/3ji7i5 \ 

;7 2 y ;7d J 



(6.8) 



Because we consider the stress-energy tensor of the point particles, we can evaluate the RHS of (6.8) explicitly, and 
we have 



rpaj ;fj _|_ rp/3j :ot }_ (bk)af3rp-fS 



= E 



dn 



;7 2 
' D 2 z? 



;7<5 



dr 



D 2 z^ 
dr? 



1 (bk)a0 D2z l yjl 

2 9 dr? V 



5W(x - Zi {n)) 



(6.9) 



where m,i, Tj and zf are the mass, the proper time and the orbital coordinates of the ith particle. Here, (D 2 zf /dr 2 ) 
is the 4-acceleration of the ith particle, and it depends on (d/dt)£, not on (d/dt)C. We consider integration of (6.8) 
using the adiabatic expansion method in §V. The result is not exactly equal to A"* 3 (AM 1 ' )), because the procedure 
of the adiabatic expansion does not commute with the derivative. However, to leading order in A, the adiabatic 
expansion and the derivative operation commute. Thus, up to 0(A 2 ), A al3 (\h^ 1,0 ^) does not contain (d/dt)C in an 
explicit manner. Because the source term of (6.7) does not have explicit dependence on the elapsed time of (d/df)C, 
we avoid the Burke problem for /i( 2 ' )"' 3 . 

Thus, we see that the source terms for h^ 1 ' 1 ^ 13 and h^ 2 - ^ 13 , the metric perturbations of 0(A 2 ) with this specific 
gauge fixing term do not depend on (d/dt)C, which grows linearly in the elapsed time after we switch on the radiation 
reaction. Although the precise behavior may still depend on the detailed form of the background metric, these source 
terms induce metric perturbations that do not grow in time, and the perturbation expansion in A holds for a sufficiently 
long time, i.e. longer than the radiation reaction time. 
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VII. SUMMARY AND DISCUSSION 



In this paper, we have investigated the so-called gauge problem arising in the study of the gravitational radiation 
reaction problem to a multi-particle system. Because the radiation reaction time scale is longer than the dynamical 
time scale of the system, we consider the use of a perturbation method to study the effect of the radiation reaction. 

Due to the radiation reaction, the evolution of the system is dissipative. There are some special perturbation 
techniques that can be applied o study a dissipative system, for example, a damped harmonic oscillator. The key 
problem that arises when using the standard perturbation expansion is that the perturbation expansion does not 
remain valid as the system evolves. Because the dissipative effect accumulates in time, the effect of the perturbation 
could become dominant over the background, and therefore the perturbation calculation may lose its predictive ability. 
In the case of a gravitational perturbation, such a problem could result from an inappropriate choice of the gauge 
condition. In Sec. Ill, we first pointed out this problem in the study of the radiation reaction given in Ref. [7], in which 
we call it the "Burke problem" . A crucial consequence of this problem is that the evolution of the system becomes 
non-adiabatic despite the fact that the prediction under the adiabatic approximation is widely accepted, because of 
the observational evidence of the Hulse- Taylor binary. 

We considered what condition must be imposed on the gauge in order to avoid this problem. For this purpose, we 
proposed a new metric perturbation scheme in §IV. The problem of the standard metric perturbation was discussed 
in Ref. [7]. Because the linearized Einstein operator is algebraically divergence free, if we simply chooose the stress- 
energy tensor as the source term of the linearized Einstein tensor, the stress-energy tensor must be conserved in the 
background, and we cannot consistently take into account the radiation reaction effect on the matter to leading order. 
By adding a gauge fixing term, we can avoid this problem. 

Here we comment that adding the gauge fixing term is not what we usually do in the black hole perturbation 
calculation, such as the Zcrilli-Rcgge- Wheeler formalism. [13] The Zerilli-Rcgge- Wheeler formalism is a convenient 
method to calculate a linear metric perturbation of a Schwarzschild black hole. In this formalism, we do not add 
a gauge fixing term to the linearized Einstein equation. Instead, we explicitly fix some components of the metric 
perturbation. In this procedure, the stress-energy tensor as a source of the metric perturbation must satisfy the 
conservation law in the background, because we only have six components of the metric perturbation. We could 
define the Green function of the Zerilli-Regge- Wheeler formalism. However, it would contain only the propagation of 
the gravitational mode, because of the explicit gauge fixing. Hence, this Green function would not be sufficient, as a 
Green function including propagation of the vector and scalar modes is necessary in this new perturbation scheme. 
The gauge fixing term in this new perturbation scheme does not vanish, and one can easily see that it comes from the 
vector and scalar modes of the Green function, because it is induced by T a/3 ;( 3. Therefore, it is necessary to extend 
the Zcrilli-Rcgge- Wheeler formalism in order to apply it to the new metric perturbation scheme. We shall discuss this 
problem elsewhere. 

Although one may consider the dissipative effect of the particles' motion to lower order in the metric perturbation 
using the new metric perturbation scheme, a systematic calculation of higher-order metric perturbations cannot be 
derived from lower-order ones. In §V, we propose another perturbation method to avoid this technical problem, which 
was originally discussed in Ref. [7] . 

Using this new metric perturbation scheme, we turned to the Burke problem with a specific choice of the gauge 
fixing term in §VI. For a complete solution of the Burke problem, it is necessary to determine the behavior of every 
order in the perturbation expansion of the metric. However, this is a difficult task. In §VI, we only computed the 
metric perturnation to second order. We further supposed that the Green function is well-behaved in the sense that, 
if the source term does not grow in time, the resulting field will not grow in time. Under this condition, we find that 
our choice of the gauge fixing term allows us to avoid the Burke problem. 
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